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On the basis of the adiabatic self-consistent collective coordinate method, we develop an effi- 
cient microscopic method of deriving the five-dimensional quadrupole collective Hamiltonian and 
illustrate its usefulness by applying it to the oblate-prolate shape coexistence/mixing phenomena in 
proton-rich 68,70 ' 72 Se. In this method, the vibrational and rotational collective masses (inertial func- 
tions) are determined by local normal modes built on constrained Hartree-Fock-Bogoliubov states. 
Numerical calculations are carried out using the pairing-plus-quadrupole Hamiltonian including the 
quadrupole-pairing interaction within the two major-shell active model spaces both for neutrons and 
protons. It is shown that the time-odd components of the moving mean-field significantly increase 
the vibrational and rotational collective masses in comparison with the Inglis-Belyaev cranking 
masses. Solving the collective Schrodinger equation, we evaluate excitation spectra, quadrupole 
transitions and moments. Results of the numerical calculation are in excellent agreement with re- 
cent experimental data and indicate that the low-lying states of these nuclei are characterized as 
an intermediate situation between the oblate-prolate shape coexistence and the so-called 7 unstable 
situation where large-amplitude triaxial-shape fluctuations play a dominant role. 

PACS numbers: 21.60.Ev, 21.10.Re, 27.50.+e 



I. INTRODUCTION 

The major purpose of this paper is to develop an effi- 
cient microscopic method of deriving the five-dimensional 
(5D) quadrupole collective Hamiltonian 1-4] and illus- 
trate its usefulness by applying it to the oblate-prolate 
shape coexistence/mixing phenomena in proton-rich Se 
isotopes ■ ^ s is well known, the quadrupole collec- 
tive Hamiltonian, also called the general Bohr-Mottelson 
Hamiltonian, contains six collective inertia masses (three 
vibrational masses and three rotational moments of in- 
ertia) as well as the collective potential. These seven 
quantities are functions of the quadrupole deformation 
variables j3 and 7, which represent the magnitude and 
triaxiality of the quadrupole deformation, respectively. 
Therefore, we also call the collective inertial masses 'in- 
ertial functions.' They are usually calculated by means 
of the adiabatic perturbation treatment of the moving 
mean field and the version taking into account nu- 
clear superfluidity is called the Inglis-Belyaev (IB) 
cranking mass or the IB inertial function. Its insuffi- 
ciency has been repeatedly emphasized, however (see e.g., 
Refs. fTlT - fl4| ). The most serious shortcoming is that the 
time-odd terms induced by the moving mean field are 
ignored, which breaks the self-consistency of the theory 
In fact, one of the most important motives of 
constructing microscopic theory of large-amplitude col- 
lective motion was to overcome such a shortcoming of 
the IB cranking mass (l5j . 

As fruits of long-term efforts, advanced microscopic 
theories of inertial functions are now available (see Refs. 
[l5H2(l ] for original papers and Refs. [13, HI] for reviews). 



These theories of large-amplitude collective motion have 
been tested for schematic solvable models and ap plied 
to heavy-ion collisions and giant resonances [l8l . |26| . For 
nuclei with pairing correlations, Dobaczewski and Skalski 
studied the quadrupole vibrational mass with use of the 
adiabatic time-dependent Hartree-Fock-Bogoliubov (AT- 
DHFB) theory and concluded that the contributions from 
the time-odd components of the moving mean-field sig- 
nificantly increase the vibrational mass compared to the 
IB cranking mass [jq |. Somewhat surprisingly, however, 
to the best of our knowledge, the ATDHFB vibrational 
masses have never been used in realistic calculations for 
low-lying quadrupole spectra of nuclei with superfluid- 
ity. For instance, in recent microscopic studies |29l - f34| by 
means of the 5D quadrupole Hamiltonian, the IB crank- 
ing formula are still used in actual numerical calcula- 
tion for vibrational masses. This situation concerning the 
treatment of the collective kinetic energies is in marked 
contrast with the remarkable progress in microscopic cal- 
culation of the collective potential using modern effective 
interactions or energy density functionals (see Ref. [35[ 
for a review). 

In this paper, on the basis of the adiabatic self- 
consistent collective coordinate (ASCC) method [36| . 
we formulate a practical method of deriving the 5D 
quadrupole collective Hamiltonian. The central con- 
cept of this approach is local normal modes built on 
constrained Hartrcc-Fock-Bogoliubov (CHFB) states [13] 
defined at every point of the (/3,7) deformation space. 
These local normal modes are determined by the lo- 
cal QRPA (LQRPA) equation that is an extension of 
the well-known quasiparticle random-phase approxima- 
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tion (QRPA) to non-equilibrium HFB states determined 
by the CHFB equations. We therefore use an abbrevia- 
tion 'CHFB+LQRPA method' for this approach. This 
method may be used in conjunction with any effec- 
tive interaction or energy density functional. In this 
paper, however, we use, for simp licity, the pairing- 
plus-quadrupole (P+Q) force [la [3!| including the 
quadrupole-pairing force. Inclusion of the quadrupole- 
pairing force is essential because it produces the time-odd 
component of the moving field [4Q| . 

To examine the feasibility of the CHFB+LQRPA 
method, we apply it to the oblate-prolate shape coex- 
istence/mixing phenomena in proton-rich 68 > 70 > 72 Se 
B EH, EH . These phenomena are taken up because we 
obviously need to go beyond the traditional framework 
of describing small-amplitude vibrations around a single 
HFB equilibrium point to describe them; that is, they 
are very suitable targets for our purpose. We shall show 
in this paper that this approach successfully describes 
large-amplitude collective vibrations extending from the 
oblate to the prolate HFB equilibrium points (and vice 
versa). In particular, it will be demonstrated that we 
can describe very well the transitional region between the 
oblate-prolate shape coexistence and the 7 unstable sit- 
uation where large-amplitude triaxial-shape fluctuations 
play a dominant role. 

This paper is organized as follows. In Sec. UH we 
formulate the CHFB+LQRPA as an approximation of 
the ASCC method and derive the 5D quadrupole col- 
lective Hamiltonian. In Sec. IIII1 we calculate the vibra- 
tional and rotational masses by solving the LQRPA equa- 
tions, and discuss their properties in comparison with 
those calculated by using the IB cranking formula. In 
Sec. |TVl we calculate excitation spectra, B(E2), and 
spectroscopic quadrupole moments of low-lying states 
in 68 > 70 > 72 Se and discuss properties of the oblate-prolate 
shape coexistence/mixing in these nuclei. Conclusions 
are given in Sec. E] 

II. MICROSCOPIC DERIVATION OF THE 5D 
QUADRUPOLE COLLECTIVE HAMILTONIAN 

A. The 5D quadrupole collective Hamiltonian 

Our aim in this section is to formulate a practical 
method of microscopically deriving the 5D quadrupole 
collective Hamiltonian [l|4l| 

W co u = T vib + T rot +V(0, 7), (1) 

T vib = \d^{P, 7)/? 2 + D M {P, 7 )/3 7 + i# 77 (/3, 7)7", 

(2) 

1 3 

T m t = ^JkWnH, (3) 
fe=i 

starting from an effective Hamiltonian for finite many- 
nucleon systems. Here, T v ;b and T ro t denote the ki- 



netic energies of vibrational and rotational motions, while 
V {(},"{) represents the collective potential. The veloci- 
ties of the vibrational motion are described in terms of 
the time-derivatives (/3, 7) of the quadrupole deforma- 
tion variables (/3, 7) representing the magnitude and the 
triaxiality of the quadrupole deformation, respectively. 
The three components Uk of the rotational angular ve- 
locity are defined with respect to the intrinsic axes asso- 
ciated with the rotating nucleus. The inertial functions 
for vibrational motions (vibrational masses), Dpp, D^, 
and -D 77 , and the rotational moments of inertia J), are 
functions of /3 and 7. 

As seen in the recent review by Prochniak and Ro- 
hozinski there are numerous papers on microscopic 
approaches to the 5D quadrupole collective Hamiltonian; 
among them, we should quote at least early papers by 
Belyaev 0, Baranger-Kumar [43|, |44|, Pomorski et al. 
0> lisf- a nd recent papers by Girod et al. [33[ , Niksic 
et al. |29ll30| . and Li et al. [HI HI]- In all these works, 
the IB cranking formula is used for the vibrational iner- 
tial functions. Below, we outline the procedure of deriv- 
ing the vibrational and rotational inertial functions on 
the basis of the ASCC method. 



B. Basic equations of the ASCC method 



To derive the 5D quadrupole collective Hamiltonian 
"Hcdi starting from a microscopic Hamiltonian H , we use 
the ASCC method plUsl]. This method enables us to de- 
termine a collective submanifold embedded in the large- 
dimensional TDHFB configuration space. We can use 
this method in conjunction with any effective interaction 
or energy density functional to microscopically derive the 
collective masses taking into account time-odd mean-field 
effects. For our present purpose, we here recapitulate 
a two-dimensional (2D) version of the ASCC method. 
We suppose existence of a set of two collective coordi- 
nates (g 1 , g 2 ) that has a one-to-one correspondence to the 
quadrupole deformation variable set (f3,-f) and try to de- 
termine a 2D collective hypersurface associated with the 
large-amplitude quadrupole shape vibrations. We thus 
assume that the TDHFB states can be written on the 
hypersurface in the following form; 



\ ( j>{q,p,^n)) = e- i ^^ N " T) \<t>{q,p,n)) 

= e -iE^ w S w e «h.P.") |0(q)) , (4) 
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with 



and 



G(g, P ,n) = Pi&iQ) + E ^ (r) e (T) (g), (5) 



i=l,2 T=n,p 



(g) = Q A (g) + <r(g) 



+ Q%{q)ola p ] , (6) 

&«(g) = Et S A (9)4 4 + e2 A *( 9 )a,a Q ]. 

(7) 



a/3 



For a gauge-invariant description of nuclei with super- 
fluidity, we need to parametrize the TDHFB state vec- 
tors, as above, not only by the collective coordinates 
Q = (q 1 ^ 2 ) an d conjugate momenta p = (pi,P2) but 
also by the gauge angles <p = Up^ n \ip^) conjugate to 
the number variables n — (n^- n ',n^) representing the 
pairing-rotational degrees of freedom (for both neutrons 
and protons). In the above equations, Q l (q) and Q^(q) 
are infinitesimal generators which are written in terms 
of the quasiparticle creation and annihilation operators 
(aT,, a a ) locally defined with respect to the moving- frame 
HFB states \4>(q)). Note that the number operators are 
defined as = iV^ T ' — subtracting the expecta- 
tion values (-/Vq"' , Nq^ ) of the neutron and proton num- 
bers at \(f>(q)). In this paper, we use units with h = 1. 

The moving- frame HFB states \<f)(q)) and the infinites- 
imal generators Q l (q) are determined as solutions of the 
moving-frame HFB equation, 



8{4>{q)\H M {q)\<t>{q)) = 0, 
and the moving-frame QRPA equations, 

S(tf,(q)\ [H M (q),Q l (q)] - - V B* k {q)P k {q) 

1 ' J 



(8) 



1 

+ 2 



L k 



\<Kq)) = o, 



(9) 



5(<f>(q)\ [H M (q)^Pi(q)}-^ C ^)Q j iQ) 



ay 



E^ (T) i^)) = o, 



(10) 



which are derived from the time-dependent variational 
principle. Here, Hm{q) is the moving- frame Hamiltonian 
given by 

H M {q) = # - E X (T \q)N^ - £ — Q\q) (11) 



k 



with 



T fe (n\ = - Rkl ( dB li , ggij _ ggij 



(12) 



). (13) 



The infinitesimal generators -P,(q) arc defined by 

d 



P i (q)\ ( p{q)}=i- J \4,( q )) 



(14) 



with 



Pi{q) ^'E^^/i " ^/3(9)a/30a], (15) 

a/3 

and determined as solutions of the moving-frame QRPA 
equations. 

The collective Hamiltonian is given as the expectation 
value of the microscopic Hamiltonian with respect to the 
TDHFB state: 

U(q,p,ri) ={(f>(q,p,n)\H\<f>(q,p,n)) 

=vto) + E \ Bij {i)PiPi + E ^ T) (q)n {T \ 

ij t 

(16) 

where 



V{q) =H(q,p,n) 



p=0,n=0 



B ij (q) 
A (T) (9) 



& 2 n 

dpidpj 

on 



p=0,n=0 



p=0,n=0 



(17) 
(18) 
(19) 



represent the collective potential, inverse of the collective 
mass, and the chemical potential, respectively. Note that 
the last term in Eq. (|10[) can be set to zero adopting the 
QRPA gauge-fixing condition, dX^/dq 1 — [4a |. 

The basic equations of the ASCC method are invari- 
ant against point transformations of the collective coordi- 
nates (q 1 , q 2 ). The B l] (q) and Cy (q) can be diagonalized 
simultaneously by a linear coordinate transformation at 
each point of q = (q 1 , q 2 ). We assume that we can intro- 
duce the collective coordinate system in which the diag- 
onal form is kept globally. Then, we can choose, without 
losing generality and for simplicity, the scale of the col- 
lective coordinates q = (g 1 ,^ 2 ) such that the vibrational 
masses become unity. Consequently, the vibrational ki- 
netic energy in the collective Hamiltonian ()16j) is written 
as 



r vib = \ E fe) 2 = ~ E (« 1 ) 2 



(20) 
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C. CHFB+LQRPA equations 

The basic equations of the ASCC method can be solved 
with an iterative procedure. This task was successfully 
carried out for extracting a one-dimensional (ID) collec- 
tive path embedded in the TDHFB configuration space 
[46l |47| . To determine a 2D hypersurface, however, 
the numerical calculation becomes too demanding at the 
present time. We therefore introduce practical approxi- 
mations as follows: First, we ignore the curvature terms 
(the third terms in Eqs. © and (JTOj) ) , which vanish at 
the HFB equilibrium points where dV/dq l = 0, assum- 
ing that their effects are numerically small. Second, we 
replace the moving- frame HFB Hamiltonian H m (<z) and 
the moving -frame HFB state \<j>{q\q 2 )) with a CHFB 
Hamiltonian -Hchfb 7) and a CHFB state \<fi(l3 , j)} , 
respectively, on the assumption that the latters are good 
approximations to the formers. 

The CHFB equations are given by 



*<M8,7)|#chfb(MI0G8,7)> =0, 



(21) 



# chfb (ft 7) =H - Y, A<T) (^l) N(T) 

T 

- Y MftT)^ (22) 



m=0,2 



with four constraints 

(^,7)1^1^/3,7)) =^ T) , (T = n,p) (23) 
W,i)\D { £\<l>(M))=D ( £>, (m = 0,2) (24) 

where -D^m denotes Hermitian quadrupole operators, 
D20 and (D22 + -C*2-2)/2 for m = and 2, respectively 
(see Ref. [46| for their explicit expressions). We define 
the quadrupole deformation variables (/3, 7) in terms of 
the expectation values of the quadrupole operators: 



/3 cos 7 = i|D^ = 77 (</>(/?, 7) I Z^ J \4>{P, 7)) , (25) 



(+) 



^=/3sin 7 = ^+ ) =^{008, 7 )|i)^|^, 7)), (26) 
where 77 is a scaling factor (to be discussed in subsection 

MM- 

The moving frame QRPA equations, (J9j) and (TlOl) . then 
reduce to 

5(^/3,7)1 [iW B (/3, 7), W, 7)] 

-ip 4 (/3, 7 )|0(/3,7)) =0, (i = l,2) (27) 
1 

and 



(+) 



<5^(/3,7)|[ff C HFB(/3,7),-A(/3,7)] 
-a(/3, 7 )Q l (/3,7)l0(/3,7)) =0. (* = 1,2) (28) 



Here the infinitesimal generators, Q 1 ^,^) and Pj(/3, 7), 
are local operators defined at (/?, 7) with respect to the 
CHFB state |^>(/3,7)). These equations are solved at 
each point of (/3, 7) to determine Q z ((3,j), Pj(/3, 7), and 
Ci(/3>7) = ujf((3,j). Note that these equations are valid 
also for regions with negative curvature (Ci(f3,j) < 0) 
where the QRPA frequency a>j(/3,7) takes an imagi- 
nary value. We call the above equations 'local QRPA 
(LQRPA) equations'. There exist more than two solu- 
tions of LQRPA equations ([27| and (|28|). and we need to 
select relevant solutions. A useful criterion for selecting 
two collective modes among many LQRPA modes will be 
given in subsection 1111 C I with numerical examples. Con- 
cerning the accuracy of the CHFB+LQRPA approxima- 
tion, some arguments will be given in subsection 1111 Fl 



D. Derivation of the vibrational masses 

Once the infinitesimal generators Q l (/3, 7) and Pi(j3, 7) 
are obtained, we can derive the vibrational masses ap- 
pearing in the 5D quadrupole collective Hamiltonian (TTJ) . 
We rewrite the vibrational kinetic energy T v ^ given by 
Eq. (|20|) in terms of the time-derivatives, and 7, of the 
quadrupole deformation variables in the following way. 
We first note that an infinitesimal displacement of the 
collective coordinates (q , q 2 ) brings about a correspond- 
ing change, 

dD ^=T,^W L ^ (™ = °>2), (29) 

i=l,2 ^ 

in the expectation values of the quadrupole operators. 
The partial derivatives can be easily evaluated as 

dD ™ {mi)\D^\mi)) 



dq % dq l 



:(0(A 7 )|[£>(+ ) ) p i (A7)]|0(/3,7)) ) (30) 



= (mi)\ J ACS, 7)] 1^08, 7)> . (31) 

without need of numerical derivatives. Accordingly, the 
vibrational kinetic energy can be written 

Tvib = \M Q0 [D { £f + M m b<£>b<£ + \m 22 [d { + ] ] 2 , 



with 



M w (/3, 7 ) = 



dq l dq % 
^ 8D (+) BD {+) 

=1,2 UU 2m alJ 2m' 



(32) 



(33) 



Taking time-derivative of the definitional equations of 
(f3, 7), Eqs. (|2"5")) and (|2^|) . we can straightforwardly trans- 
form the above expression (|32[) to the form in terms of 
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(0, 7). The vibrational masses (Dpp, Dp 7 , -D77) are then 
obtained from (Moo, M 02 , M 22 ) through the following re- 
lations: 



D/30=r] 2 [Moo cos 2 7 + V2M Q2 sin 7 cos ' 



+ -M 22 sm z 7j , 

Dfs-y =j3-q- 2 [-M00 sin 7 cos 7 



(34) 



— ^=Af 2(cos 2 7 — sin 2 7) + \-M22 sin 7 cos 7 
V 2 2 



£> 77 =/3 2 ry 2 ^M o sin 2 7 - V2Af 2 sin 7 cos 7 



+ 2 M 22 cos 2 7 



(35) 



(36) 



E. Calculation of the rotational moments of inertia 

We calculate the rotational moments of inertia Jk (0, 7) 
using the LQRPA equation for the collective rotation [46[ 
at each CHFB state, 



*W,7)IPc 



HFB, ^fc 



1 



T (Jk)- 1 I k \cl>W,7))=0, (37) 



(^C8,7)| [**(/9, 7),^] 1008, 7)> = (38) 

where ^fk(0, 7) and represent the rotational angle and 
the angular momentum operators with respect to the 
principal axes associated with the CHFB state \<f>(0, 7)). 
This is an extension of the Thouless-Valatin equation J48[ 
for the HFB equilibrium state to non-equilibrium CHFB 
states. The three moments of inertia can be written as 

Jit 08, 7) = 40 2 D k (0, 7) sin 2 7fc (fc = 1,2, 3) (39) 

with 7^ = 7 — (27rfe/3). If the inertial functions Dk(0, 7) 
above are replaced with a constant, then 7) re- 

duce to the well-known irrotational moments of inertia. 
In fact, however, we shall see that their (/3, 7) depen- 
dence is very important. We call Jk{0il) and 7) 
determined by the above equation 'LQRPA moments of 
inertia' and 'LQRPA rotational masses', respectively. 



F. Collective Schrodinger equation 

Quantizing the collective Hamiltonian (|TJ) with the 
Pauli prescription, we obtain the collective Schrodinger 
equation Q 

{Tvib + Trot + V}* aIM {0, 7, 0) = E a i* aIM (0, 7, si), 

(40) 



where 



1 



1 



2VWR 4 



d 2 



R 



1 



2 sin 37 



3 i 2 



sin 3^/Dp 1 dj3 



— SW ->>~D ■„.;(). 



(41) 
(42) 



with 

i?(/3, 7 ) =D 1 (0, 1 )D 2 {0, 1 )D z {0, 1 ), 
W(0,j)={D^(0, 7 )D^(0 n ) - [D 01 {0, 1 )] 2 }0 



(43) 

-2 

(44) 



The collective wave function in the laboratory frame, 
^aiM(0,l,to), is a function of 0, 7, and a set of three 
Euler angles fi. It is specified by the total angular mo- 
mentum /, its projection onto the z-axis in the laboratory 
frame M, and a that distinguishes the eigenstates pos- 
sessing the same values of / and M. With the rotational 
wave function D I IiIK (n), it is written as 



aIK(0,j)(to\IMK), (45) 

K— even 

where 

(to\iMK) = J—^±L—p^ K (Q) + (-yv^m 



16tt 2 (1 + 4o)' 



(46) 



The vibrational wave functions in the body-fixed frame, 
^aiK{0, 7), are normalized as 



d0d 1 \^ aI (0, 1 )\ 2 \G(0, 7 )\- 



1. 



where 



|$ q7 (/3, 7 )| 2 ^ \**iK(fi,i)? 



(47) 



(48) 



and the volume element \G(0,j)\^d0d / y is given by 

\G{0, 1 )\?d0d 1 = 20 i ^/W(0, 7 )i?(/3, 7) sin 3 7 d/3d 7 . 

(49) 

Thorough discussions of their symmetries and the bound- 
ary conditions for solving the collective Schrodinger equa- 
tion are given in Refs. [l|-0|- 
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III. CALCULATION OF THE COLLECTIVE 
POTENTIAL AND THE COLLECTIVE MASSES 



A. Details of numerical calculation 



The CHFB+LQRPA method outlined in the preced- 
ing section may be used in conjunction with any effec- 
tive interaction, e.g., density-dependent effective inter- 
actions like Skyrme forces, or modern nuclear density 
functionals. In this paper, as a first step toward such 
calculations, we use a version of the P+Q force model 
38, 39] that includes the quadrupole-pairing in addi- 
tion to the monopole-pairing interaction. Inclusion of 
the quadrupole-pairing is essential, because neither the 
monopole-pairing nor the quadrupole particle-hole inter- 
action contributes to the time-odd mean-field effects on 
the collective masses (lfj; that is, only the quadrupole- 
pairing induces the time-odd contribution in the present 
model. Note that the quadrupole-pairing effects were 
not considered in Ref. |l6j. In the numerical calculation 
for 68 > 70 > 72 Se presented below, we use the same notations 
and parameters as in our previous work [47]. The shell 
model space consists of two major shells (iV s h = 3, 4) for 
neutrons and protons and the spherical single-particle en- 
ergies are calculated using the modified oscillator poten- 
tial [i^, [5(| • The monopole-pairing interaction strengths 

(for neutrons and protons), and the quadrupole- 

particle-hole interaction strength, are determined such 
that the magnitudes of the quadrupole deformation /3 and 
the monopole-pairing gaps (for neutrons and protons) 
at the oblate and prolate local minima in 68 Se approxi- 
mately reproduce those obtained in the Skyrme-HFB cal- 
culations [Hlj]. The interaction strengths for 70 Se and 
72 Se are then determined assuming simple mass-number 
dependence G ( q t) ~ A -1 and X ' = ~ A~i 

(b denotes the oscillator- length parameter). For the 
quadrupole-pairing interaction strengths (for neutrons 
and protons), we use the Sakamoto-Kishimoto prescrip- 
tion [52j to derive the self-consistent values. Following 
the conventional treatment of the P+Q model (53[, we 
ignore the Fock term, so that we use the abbreviation 
HB (Hartree-Bogoliubov) in place of HFB in the follow- 
ing. In the case of the conventional P+Q model, the HB 
equation reduces to a simple Nilsson + BCS equation 
(see, e.g., Ref. H3). The presence of the quadrupole- 
pairing interaction in our case does not allow such a re- 
duction, however, and we directly solve the HB equation. 
In the P+Q model, the scaling factor rj in Eqs. (|25j) and 
([2T5)) is given by 77 = x'/^ob 2 , where ujo denotes the 
frequency of the harmonic-oscillator potential. Effective 
charges, (e„,e p ) = (0.4, 1.4), are used in the calculation 
of quadrupole transitions and moments. 

To solve the CHB + LQRPA equations on the (j3, 7) 
plane, we employ a two-dimensional mesh consisting of 
3600 points in the region < /3 < 0.6 and 0° < 7 < 60°. 



Each mesh point (/3j,7j) is represented as 

Pi = (i - 0.5) X 0.01, (i = l,---60), (50) 
7i = (j-0.5)xl°, (7 = 1,... 60). (51) 

One of the advantages of the present approach is that we 
can solve the CHB + LQRPA equations independently 
at each mesh point on the (/3, 7) plane, so that it is suited 
to parallel computation. 

Finally, we summarize the most important differences 
between the present approach and the Baranger-Kumar 
approach (43[. First, as repeatedly emphasized, we intro- 
duce the LQRPA collective massess in place of the crank- 
ing masses. Second, we take into account the quadrupole- 
pairing force (in addition to the monopole-pairing force), 
which brings about the time-odd effects on the collective 
masses. Third, we exactly solve the CHB self-consistent 
problem, Eq. (|2ip. at every point on the (/3, 7) plane us- 
ing the gradient method, while in the Baranger-Kumar 
works the CHB Hamiltonian is replaced with a Nilsson- 
like single- particle model Hamiltonian. Fourth, we do 
not introduce the so-called core contributions to the col- 
lective masses, although we use the effective charges to 
renormalize the core polarization effects (outside of the 
model space consisting of two major shells) into the 
quadrupole operators, We shall see that we can well re- 
produce the major characteristics of the experimental 
data without introducing such core contributions to the 
collective masses. Fifth, most importantly, the theoret- 
ical framework developed in this paper is quite general, 
that is, it can be used in conjunction with modern density 
functionals going far beyond the P+Q force model. 



B. Collective potentials and pairing gaps 

We show in Fig.[T]the collective potentials V((3, 7) cal- 
culated for 68,70 ' 72 Se. It is seen that two local minima 
always appear both at the oblate (7 = 60°) and prolate 
(7 = 0°) shapes, and, in all these nuclei, the oblate mini- 
mum is lower than the prolate minimum. The energy dif- 
ference between them is, however, only several hundred 
keV and the potential barrier is low in the direction of 
triaxial shape (with respect to 7) indicating 7-soft char- 
acter of these nuclei. In Fig. [T] we also show the collective 
paths (connecting the oblate and prolate minima) deter- 
mined by using the ID version of the ASCC method [47$ . 
It is seen that they always run through the triaxial valley 
and never go through the spherical shape. 

In Fig. the monopole- and quadrupole-pairing gaps 
calculated for 68 Se are displayed. They show a signifi- 
cant (/3, 7) dependence. Broadly speaking, the monopole 
pairing decreases while the quadrupole pairing increases 
as f3 increases. 
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FIG. 1: (Color online) Collective potential V(/3,j) for 
68,70,72 g e re gj ons hig ner than 3 MeV (measured from 

the oblate HB minima) are drawn by rosybrown color. One- 
dimensional collective paths connecting the oblate and prolate 
local minima are determined by using the ASCC method and 
depicted with bold red lines. 



C. Properties of the LQRPA modes 

In Fig. |3] the frequencies squared, wf(/3, 7), of various 
LQRPA modes calculated for 68 Se are plotted as func- 
tions of /? and 7. In the region of the (j3, 7) plane where 
the collective potential energy is less than about 5 MeV, 
we can easily identify two collective modes among many 
LQRPA modes, whose ujf((3,^f) are much lower than 
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FIG. 2: (Color online) Monopole- and quadrupole-pairing 
gaps for neutrons of 68 Se are plotted in the (/?, 7) deforma- 
tion plane. ( upper left) Monopole pairing gap Aq \ (lower 
left) Quadrupole pairing gap . (lower right) Quadrupole 



pairing gap A^™' • 



See Rcf. 



for definitions of Ai n \ A™ , 



and A 



(n) 



those of other modes. Therefore we adopt the two low- 
est frequency modes to derive the collective Hamiltonian. 
This result of numerical calculation supports our assump- 
tion that there exists a 2D hypersurface associated with 
large-amplitude quadrupole shape vibrations, which is 
approximately decoupled from other degrees of freedom. 
The situation changes when the collective potential en- 
ergy exceeds about 5 MeV and/or the monopole-pairing 
gap becomes small. A typical example is presented in the 
bottom panel of Fig. [3] It becomes hard to identify two 
collective modes well-separated from other modes when 
j3 > 0.4, where the collective potential energy is high (see 
Fig. [IJ and the monopole-pairing gap becomes small (see 
Fig. [2]) . In this example, the second- lowest LQRPA mode 
in the 0.4 < /3 < 0.5 region has pairing-vibrational char- 
acter but becomes non-collective for j3 > 0.5. In fact, 
many non-collective two-quasiparticle modes appear in 
its neighborhood. This region in the (/3, 7) plane is not 
important, however, because only tails of the collective 
wave function enter into this region. 

It may be useful to set up a prescription that works 
even in a difficult situation where it is not apparent how 
to choose two collective LQRPA modes. We find that the 
following prescription always works well for selecting two 
collective modes among many LQRPA modes. This may 
be called a minimal metric criterion. At each point on 
the (/3,7) plane, we evaluate the vibrational part of the 
metric W(/3, 7) given by Eq. (|44p for all combinations of 
two LQRPA modes, and find the pair that gives the mini- 
mum value. We show in Fig.Uhow this prescription actu- 
ally works. In this figure, the W(/3, 7) values are plotted 
as functions of j3 and 7 for many pairs of the LQRPA 
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FIG. 3: Frequencies squared lo 2 of the LQRPA modes calcu- 
lated for 68 Se are plotted as functions of /3 or 7. The LQRPA 
modes adopted for calculation of the vibrational masses are 
connected with solid lines, (top) Dependence on 7 at /3 = 0.3. 
(middle) Dependence on j3 along the 7 = 0.5° line, (bottom) 
Dependence on /3 along the 7 = 30.5° line. 



modes. In the situations where the two lowest-frequency 
LQRPA modes are well separated from other modes, this 
prescription gives the same results as choosing the two 
lowest-frequency modes (see the top and middle panels). 
On the other hand, a pair of the LQRPA modes different 
from the lowest two modes is chosen by this prescrip- 
tion in the region mentioned above (the bottom panel). 
This choice may be better than that using the lowest- 
frequency criterion, because we often find that a nor- 
mal mode of pairing vibrational character becomes the 
second lowest LQRPA mode when the monopole-pairing 
gap significantly decreases in the region of large /?. The 
small values of the vibrational metric implies that the di- 



FIG. 4: Dependence on /? and 7 of the vibrational part of the 
metric W(/3,y) calculated for 68 Se. (top) Dependence on 7 at 
/? = 0.3. (mtddle) Dependence on /3 along the 7 = 0.5° line. 
(bottom) Dependence on /3 along the 7 = 30.5° line. The cross 
symbols indicate values of the vibrational metric calculated 
for various choices of two LQRPA modes from among the 
lowest 40 LQRPA modes; the lowest mode is always chosen 
and the other is from the remaining 39 modes. The smallest 
vibrational metric is shown by solid line. For reference, the 
vibrational metric calculated using the IB vibrational mass is 
indicated by broken lines. 

rection of the infinitesimal displacement associated with 
the pair of the LQRPA modes has a large projection onto 
the (/3, 7) plane. Therefore, this prescription may be well 
suited to our purpose of deriving the collective Hamilto- 
nian for the (/3, 7) variables. It remains as an interesting 
open question for future to examine whether or not the 
explicit inclusion of the pairing vibrational degree of free- 
dom as another collective variable will give us a better 
description in such situations. 
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FIG. 5: (Color online) Vibrational masses, Dpp{p,"i), 
Dfr(B,i)/j3, and D 77 (/3, 7 )//3 2 , in unit of MeV" 1 calculated 
for 68 Se. 



D. Vibrational masses 

In Fig. [5] the vibrational masses calculated for 68 Se are 
displayed. We see that their values exhibit a significant 
variation in the (/?, 7) plane. In particular, the increase 
in the large /3 region is remarkable. 

Figure [S] shows how the ratios of the LQRPA vibra- 
tional masses to the IB vibrational masses vary on the 
(/3, 7) plane. It is clearly seen that the LQRPA vibra- 
tional masses are considerably larger than the IB vibra- 
tional masses and their ratios change depending on (3 
and 7. In this calculation, the IB vibrational masses are 




0.1 0.2 0.3 0.4 0.5 0.6 

P 

FIG. 6: (Color online) Ratios of the LQRPA vibra- 
tional masses to the IB vibrational masses, Dpp/D^ and 
D 77 /D^ I f' , calculated for 68 Se. 



evaluated using the well-known formula: 



(H 



dH, 



CHB 



8H ( 



CHB 



[£„(/?, 7) 7)F 



(52) 
(6 = P or 7 ) 

where £^(,8,7), |0), and \jjLV) denote the quasiparticle 
energy, the CHB state \(j>(/3, 7)) and the two-quasiparticle 
state atdp \(j>(fi, 7)), respectively (see Ref. [HI for the 
meaning of the indices /i and v). 

The vibrational masses calculated for 70 ' 72 Se exhibit 
behaviors similar to those for 68 Se. 



E. Rotational masses 

In Fig. [3 the LQRPA rotational masses Dk(P,j) cal- 
culated for 68 Se are displayed. Similarly to the vibra- 
tional masses discussed above, the LQRPA rotational 
masses also exhibit a remarkable variation over the (/?, 7) 
plane, indicating a significant deviation from the irrota- 
tional property. 
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FIG. 7: (Color online) Rotational masses _Dfc(/3,7) in unit of 
MeV" 1 , calculated for 68 Se. See Eq. (39]) for the relation with 
the rotational moments of inertia Jk{fi, y). 




0.1 0.2 0.3 0.4 0.5 0.6 

P 




0.1 0.2 0.3 0.4 0.5 0.6 

P 

FIG. 8: (Color online) Ratios of the LQRPA rotational masses 
to the IB rotational masses, D k (P,y)/D ( ™ ) (P,y), calculated 



Figure [8] shows how the ratios of the LQRPA rotational 

(IB) 

masses 7) to the IB cranking masses D k (/?, 7) 

vary on the (/?, 7) plane. The rotational masses calcu- 
lated for 70 ' 72 Se exhibit behaviors similar to those for 
68 Se. 

As we have seen in Figs. [5HS1 not only the vibra- 
tional and rotational masses but also their ratios to the 
IB cranking masses exhibit an intricate dependence on 
P and 7. For instance, it is clearly seen that the ra- 

(IB) 

tios, Dk(P,^)/D y k (/3,7), gradually increase as p de- 
creases. This result is consistent with the calcula- 
tion by Hamamoto and Nazarewicz (54|, where it is 



shown that the ratio of the Migdal term to the crank- 
ing term in the rotational moment of inertia (about the 
1st axis) increases as /3 decreases. Needless to say, the 
Migdal term (also called the Thouless-Valation correc- 
tion) corresponds to the time-odd mean-field contribu- 
tion taken into account in the LQRPA rotational masses, 
so that the result of Ref. implies that the ratio 

(IB) 

Di(P,^)/D\ '(/3,y), increases as p decreases, in agree- 
ment with our result. To understand this behavior, it is 
important to note that, in the present calculation, the dy- 
namical effect of the time-odd mean-field on £>i(/3,y) is 
associated with the K = 1 component of the quadrupole- 
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pairing interaction and it always works and increase the 
rotational masses, in contrast to the behavior of the static 
quantities like the magnitude of the quadrupole-pairing 
gaps, A20 and A22, which diminish in the spherical shape 
limit. Obviously, this qualitative feature holds true ir- 
respective of details of our choice of the monopole and 
quadrupole pairing interaction strengths. 

The above results of calculation obviously indicate the 
need to take into account the time-odd contributions 
to the vibrational and rotational masses by going be- 
yond the IB cranking approximation. In Refs. }29l - l32j . 
a phenomenological prescription is adopted to remedy 
the shortcoming of the IB cranking masses; that is, a 
constant factor in the range 1.40-1.45 is multiplied to 
the IB rotational masses. This prescription is, however, 
insufficient in the following points. First, the scaling 
only of the rotational masses (leaving the vibrational 
masses aside) violates the symmetry requirement for the 
5D collective quadrupole Hamiltonian [lH3| (a similar 
comment is made in Ref. [3]). Second, the ratios take 
different values for different LQRPA collective masses 
(Dpp, -D/3 7 , -D-y-y, D\, D2, and D3). Third, for every col- 
lective mass, the ratio exhibits an intricate dependence 
on /3 and 7. Thus, it may be quite insufficient to sim- 
ulate the time-odd mean-field contributions to the col- 
lective masses by scaling the IB cranking masses with a 
common multiplicative factor. 



F. Check of self-consistency along the collective 
path 

As discussed in SeaQll the CHB+LQRPA method is a 
practical approximation to the ASCC method. It is cer- 
tainly desirable to examine the accuracy of this approxi- 
mation by carrying out a fully self-consistent calculation. 
Although, at the present time, such a calculation is too 
demanding to carry out for a whole region of the (/3, 7) 
plane, we can check the accuracy at least along the ID 
collective path. This is because the ID collective path is 
determined by carrying out a fully self-consistent ASCC 
calculation for a single set of collective coordinate and 
momentum. The ID collective paths projected onto the 
(/?, 7) plane are displayed in Fig. [T] Let us use a nota- 
tion \<fc(q)) for the moving- frame HB state obtained by 
self-consistently solving the ASCC equations for a single 
collective coordinate q [HI, HtJ. To distinguish from it, 
we write the CHB state as \4>(/3(q), 7(9))). This notation 
means that the values of (3 and 7 are specified by the col- 
lective coordinate q along the collective path. In other 
words, \(f)((3(q), 7(g))) has the same expectation values of 
the quadrupole operator as those of \(f>{q)) ■ It is impor- 
tant to note, however, that they are different from each 
other, because |</>(/3(g), 7(g))) is a solution of the CHB 
equation which is an approximation of the moving-frame 
HB equation. Let us evaluate various physical quantities 
using the two state vectors and compare the results. 

In Fig. [S] various physical quantities (the pairing gaps, 
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FIG. 9: (Color online) Comparison of physical quantities eval- 
uated with the CHB + LQRPA approximation and those with 
the ASCC method. Both calculations are carried out along 
the ID collective path for 68 Se and the results are plotted as 
function of 7(g). From the top to the bottom: 1) the col- 
lective potential, 2) monopole-pairing gaps, Aq"" and Ag P \ 
for neutrons and protons, 3) frequencies squared uj 2 of the 
lowest and the second-lowest modes obtained by solving the 
moving-frame QRPA and the LQRPA equations, 4) vibra- 
tional masses, Dpp, D/s^//3, and D 11 j '/3 2 , 5) rotational masses 
Dk- In almost all cases, results of the two calculations are in- 
distinguishable, because they agree within the widths of the 
line. 



the collective potential, the frequencies of the local 
normal modes, the rotational masses, and vibrational 
masses) calculated using the moving-frame HB state 
\(j){q)} and the CHB state |0(/3(g) , 7(g))) are presented 
and compared. These calculations are carried out along 
the ID collective path for 68 Se. Apparently, the results 
of the two calculations are indistinguishable in almost all 
cases, because they agree within the widths of the line. 
This good agreement implies that the CHB+LQRPA is 
an excellent approximation to the ASCC method along 
the collective path on the (/3, 7) plane. As we shall see 
in the next section, collective wave functions distribute 
around the collective path. Therefore, it may be reason- 
able to expect that the CHB+LQRPA method is a good 
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approximation to the ASCC method and suited, at least, 
for describing the oblate-prolate shape mixing dynamics 
in 68 > 70 ' 72 Se. 



IV. LARGE-AMPLITUDE SHAPE-MIXING 
PROPERTIES OF 68 ' 70 > 72 SE 

We have calculated collective wave functions solving 
the collective Schrodinger equation (|40| and evaluated 
excitation spectra, quadrupole transition probabilities, 
and spectroscopic quadrupole moments. The results for 
low-lying states in 68 > 70 > 72 Se are presented in Figs. fTTTHTSI 

In Figs. HUl [T2l and[JJJ excitation spectra and B(E2) 
values for 68 Se, 70 Se, and 72 Se, calculated with the 
CHB+LQRPA method, are displayed together with ex- 
perimental data. The eigenstates are labeled with T* = 
+ ,2+,4 + , and 6 + . In these figures, results obtained us- 
ing the IB cranking masses are also shown for the sake 
of comparison. Furthermore, the results calculated with 
the (1+3)D version of the ASCC method reported in our 
previous paper 47[ are shown also for comparison with 
the 5D calculations. We use the abbreviation (1+3)D 
to indicate that a single collective coordinate along the 
collective path describing large-amplitude vibration and 
three rotational angles associated with the rotational mo- 
tion are taken into account in these calculations. The 
classification of the calculated low-lying states into fami- 
lies of two or three rotational bands is made according to 
the properties of their vibrational wave functions. These 
vibrational wave functions are displayed in Figs. [TT] [T3J 
and Q~5] In these figures, only the /3 4 factor in the vol- 
ume element (|49p are multiplied to the vibrational wave 
functions squared leaving the sin 37 factor aside. This 
is because all vibrational wave functions look like triax- 
ial and the probability at the oblate and prolate shapes 
vanish if the sin 37 factor is multiplied by them. 

Let us first summarize the results of the CHB+LQRPA 
calculation. The most conspicuous feature of the low- 
lying states in these proton-rich Se isotopes is the dom- 
inance of the large-amplitude vibrational motion in the 
triaxial shape degree of freedom. In general, the vibra- 
tional wave function extends over the triaxial region be- 
tween the oblate (7 = 60°) and the prolate (7 = 0°) 
shapes. In particular, this is the case for the + states 
causing their peculiar behaviors; for instance, we obtain 
two excited + states located slightly below or above the 
2^~ state. Relative positions between these excited states 
are quite sensitive to the interplay of large-amplitude 7- 
vibrational modes and the /3-vibrational modes. This 
result of calculation is consistent with the available ex- 
perimental data where the excited + state has not yet 
been found, but more experimental data are needed to 
examine the validity of the theoretical prediction. Below, 
let us examine characteristic features of the theoretical 
spectra more closely for individual nuclei. 

For 68 Se, we obtain the third band in low energy. The 



O2 and 2% states belonging to this band are also shown in 
Fig. 1101 Their vibrational wave functions exhibit nodes 
in the (3 direction (see Fig. ITTj) indicating that a (3- 
vibrational mode is excited on top of the large-amplitude 
7 vibrations. As a matter of course, this kind of state is 
outside of the scope of the (1+3)D calculation. The vi- 
brational wave functions of the yrast 2\ and states 
exhibit localization in a region around the oblate shape, 
while the yrare 2%, 4 \, and 6 J states localize around the 
prolate shape. It is apparent, however, that all the wave 
functions significantly extend from 7 = 0° to 60° over 
the triaxial region, indicating 7-soft character of these 
states. In particular, the yrare 4^~ and 6^" wave functions 
exhibit two-peak structure consisting of the prolate and 
oblate peaks. The peaks of the vibrational wave function 
gradually shift toward a region of larger f3 as the angu- 
lar momentum increases. This is a centrifugal effect de- 
creasing the rotational energy by increasing the moment 
of inertia. In the (1+3)D calculation, this effect is absent 
because the collective path is fixed at the ground state. 
Thus, the 5D calculation yields, for example, a much 
larger value for B(E2; 6^ — > 4^) in comparison with the 
(1+3)D calculation. Actually, in the 5D CHB+LQRPA 
calculation, the wave function of the yrast 6f state lo- 
calizes in the triaxial region (see Fig. ITTj) where the mo- 
ment of inertia takes a maximum value. This leads to 
a small value for the spectroscopic quadrupole moment 
(see Fig.lTTJf because of the cancellation between the con- 
tributions from the oblate-like and prolate-like regions. 
This cancellation mechanism due to the large-amplitude 
7 fluctuation is effective also in other states; although the 
spectroscopic quadrupole moments of the yrast 2^ and 
4^ (yrare 2^,4^", and 6^~) states are positive (negative) 
indicating their oblate-like (prolate-like) character, their 
absolute magnitudes are rather small. 

The .E2-transition probabilities exhibit a pattern 
reminiscent of the 7-unstable situation; for instance, 
B{E2; &+ -> 6+), B{E2; 4+ 4+), and B{E2- 2+ -> 2+) 
are much larger than B(E2; 6+ -> 4+), B(E2; 4+ -> 2+), 
and B(E2;2+ -> 0+); see Fig. [TOJ Thus, the low-lying 
states in 68 Se may be characterized as an intermedi- 
ate situation between the oblate-prolate shape coexis- 
tence and the Wilets-Jean 7-unstable model [5a]. Using 
the phenomenological Bohr-Mottelson collective Hamil- 
tonian, we have shown in Ref. (56j that it is possible to 
describe the oblate-prolate shape coexistence and the 7- 
unstable situation in a unified way varying a few param- 
eters controlling the degree of oblate-prolate asymmetry 
in the collective potential and the collective masses. The 
two-peak structure seen in the 4^~ and states may 
be considered as one of the characteristics of the inter- 
mediate situation. It thus appears that the excitation 
spectrum for 68 Se (Fig. [TTJI) serves as a typical example 
of the transitional phenomena from the 7-unstable to the 
oblate-prolate shape coexistence situations. 

Let us make a comparison between the spectra in 
Fig. [10] obtained with the LQRPA collective masses and 
that with the IB cranking masses. It is obvious that 
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FIG. 10: Excitation spectra and B(E2) values calculated for 68 Se by means of the CHB+LQRPA method (denoted 
CHB+LQRPA) and experimental data [5|-[7j]. For comparison, results calculated using the IB cranking masses (denoted 
CHB+IB) and those obtained using the (1+3)D version of the ASCC method (denoted (1+3)D ASCC) are also shown. Only 
B(E2) , s larger than 1 Weisskopf unit (in the (1+3)D ASCC and/or the CHB+LQRPA calculations) are shown in units of 
e 2 fm 4 . 




FIG. 11: (Color online) Vibrational wave functions squared l3 4 \&ik(P, 7)| 2 calculated for 68 Se. 



the excitation energies are appreciably overestimated in 
the latter. This result is as expected from the too low 
values of the IB cranking masses. The result of our cal- 
culation is in qualitative agreement with the HFB-based 
configuration-mixing calculation reported by Ljungvall et 
al. [8[ in that both calculations indicate the oblate (pro- 
late) dominance for the yrast (yrare) band in 68 Se. Quite 
recently, the B(E2;2^ — > Of) value has been measured 



in experiment Q. The calculated value (492 e 2 fm 4 ) is in 
fair agreement with the experimental data (432 e 2 fm 4 ). 

The result of calculation for 70 Se (Figs. fl2l and fl3f is 
similar to that for 68 Se. The vibrational wave functions of 
the yrast 2^~, 4f, and 6^ states localize in a region around 
the oblate shape, exhibiting, at the same time, long tails 
in the triaxial direction. We note here that, differently 
from the 68 Se case, the wave function keeps the oblate- 
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(1+3)DASCC CHB+IB CHB+LQRPA EXP 

FIG. 12: Same as Fig. [I0]but for 70 Se. Experimental data is taken from Refs. lilj. 
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FIG. 13: (Color online) Same as Fig. [TT] but for 70 Se. 



like structure. On the other hand, the yrare 2^,42, and 
d>2 states localize around the prolate shape, exhibiting, at 
the same time, small secondary bumps around the oblate 
shape. For the yrare 2^ state, we obtain a strong oblate- 
prolate shape mixing in the (1+3)D calculation [47|. This 
mixing becomes weaker in the present 5D calculation, re- 
sulting in the reduction of the B(E2;A^ — > 2%) value. 
Similarly to 68 Se, we obtain two excited + states in low 
energy. We see considerable oblate-prolate shape mix- 
ings in their vibrational wave functions, but, somewhat 
differently from those in 68 Se, the second and third + 



states in 70 Se exhibit clear peaks at the oblate and prolate 
shapes, respectively, Their energy ordering is quite sen- 
sitive to the interplay of the large-amplitude 7 vibration 
and the /3 vibrational modes. The calculated spectrum 
for 70 Se is in fair agreement with the recent experimental 
data [4l[ , although the B(E2) values between the yrast 
states are overestimated. 

The result of calculation for 72 Se (Figs. [14] and [15]) 
presents a feature somewhat different from those for 68 Se 
and 70 Se; that is, the yrast 2^,4^, and 6* states local- 
ize around the prolate shape instead of the oblate shape. 
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(1+3)DASCC CHB+IB CHB+LQRPA EXP 

FIG. 14: Same as Fig. [I0]but for 72 Se. Experimental data is taken from Refs. @, El]. 
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FIG. 15: (Color online) Same as Fig. [TT] but for 72 Se. 



The localization develops with increasing angular mo- 
mentum. On the other hand, similarly to the 68 ' 70 Se 
cases, the yrare 2^~,4^, and 6 J states exhibit the two- 
peak structure. The spectroscopic quadrupole moments 
of the 2^, 4^~, and states are negative, and their abso- 
lute magnitude increases with increasing angular momen- 
tum (see Fig. fT75)) reflecting the developing prolate char- 
acter in the yrast band, while those of the yrare states 
are small because of the two-peak structure of their vi- 
brational wave functions, that is, due to the cancellation 
of the contributions from the prolate-like and oblate-like 



regions. Also for 72 Se, we obtain two excited + states 
in low energy, but they show features somewhat differ- 
ent from the corresponding excited + states in 68 ' 70 Se. 
Specifically, the vibrational wave functions of the sec- 
ond and third + states exhibit peaks at the prolate and 
oblate shape, respectively. As seen in Fig. [TH our re- 
sults of calculation for the excitation energies and B(E2) 
values are in good agreement with the recent experimen- 
tal data [1] for the yrast 2^~,4f , and 6f states in 72 Se. 
Experimental _E2-transition data are awaited for under- 
standing the nature of the observed excited band. 
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FIG. 16: Spectroscopic quadrupole moments for 68 ' 70 > 7 Se. 
Values calculated with the LQRPA collective masses are 
shown with the triangles. For comparison, values calculated 
with the IB collective masses and those obtained with the 
(1+3)D version of the ASCC method are also shown with the 
squares and the circles, respectively. The filled symbols show 
the values for the yrast states, while the open symbols those 
for the yrare states. 

V. CONCLUSIONS 

On the basis of the ASCC method, we have developed a 
practical microscopic approach, called CHFB+LQRPA, 
of deriving the 5D quadrupole collective Hamiltonian 
and confirmed its efficiency by applying it to the oblate- 
prolate shape coexistence / mixing phenomena in proton- 
rich 68,70 ' 72 Se. The results of numerical calculation for 
the excitation energies and B(E2) values are in good 
agreement with the recent experimental data @,IH for the 



yrast 2+, 4^, and 6^ states in these nuclei. It is shown 
that the time-odd components of the moving mean-field 
significantly increase the vibrational and rotational col- 
lective masses and make the theoretical spectra in much 
better agreement with the experimental data than cal- 
culations using the IB cranking masses. Our analysis 
clearly indicates that low-lying states in these nuclei pos- 
sess a transitional character between the oblate-prolate 
shape coexistence and the so-called 7 unstable situation 
where large-amplitude triaxial-shape fluctuations play a 
dominant role. 

Finally, we would like to list a few issues for the fu- 
ture that seems particularly interesting. First, fully self- 
consistent solution of the ASCC equations for determin- 
ing the two-dimensional collective hypersurface and ex- 
amination of the validity of the approximations adopted 
in this paper in the derivation of the CHFB+LQRPA 
scheme. Second, application to various kind of collec- 
tive spectra associated with large-amplitude collective 
motions near the yrast lines (as listed in Ref. |28|). 
Third, possible extension of the quadrupole collective 
Hamiltonian by explicitly treating the pairing vibrational 
degrees of freedom as additional collective coordinates. 
Fourth, use of the Skyrme energy functionals + density- 
dependent contact pairing interaction in place of the 
P+Q force, and then modern density functionals cur- 
rently under active development. Fifth, application of 
the CHFB+LQRPA scheme to fission dynamics. The 
LQRPA approach enables us to evaluate, without need 
of numerical derivatives, the collective inertia masses in- 
cluding the time-odd mean-field effects. 
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